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The energy of Skyrmions is calculated with the help of a technique based on the excitonic repre- 
sentation: the basic set of one-exciton states is used for the perturbation-theory formalism instead 
of the basic set of one-particle states. We use the approach, at which a skyrmion-type excitation 
(at zero Lande factor) is considered as a smooth non-uniform rotation in the 3D spin space. The 
result within the framework of an excitonically diagonalized part of the Coulomb Hamiltonian can 
be obtained by any ratio r c = (e 2 /£Ib)/TUo c [where e 2 /els is the typical Coulomb energy (Ib being 
the magnetic length); uj c is the cyclotron frequency], and the Landau-level mixing is thereby taken 
into account. In parallel with this, the result is also found exactly, to second order in terms of the 
rc (if supposing r c to be small) with use of the total Hamiltonian. When extrapolated to the region 
r c ~ 1, our calculations show that the skyrmion gap becomes substantially reduced in comparison 
with the Hartree-Fock calculations. This fact brings the theory essentially closer to the available 
experimental data. 

PACS numbers: 73.43. Cd, 71.27.+a 



I. INTRODUCTION 

Up to now, in two-dimensional (2D) electron gas 
(EG) the quantum Hall effect (QHE) at the Landau- 
level (LL) filling factor of v = Ijrtpi attracted much 
theoreticalETEl and experimentalcTEil attention. The 
interest is explained by the fact that in this regime 
at vanishing (or considerably reduced) Zeeman cou- 
pling peculiar Fermionic excitations exist: namely, 
these are skyrmions, characterized by a large spin, 
\SS Z \ 3> 1, and by a topological invariant (or topo- 
logical charge) in terms..©f a field theory of the clas- 
sical 2D ferromagnet .113113 The first mapping of the 
spin-polarized quantum Hall system to an appropri- 
ate nonlinear 0(3) model, of the 2D fexromagnet, 
was used in the work of Sondhi et al a to calcu- 
late the skyrmion energy. Within the context of 
the phenomenological approach employed the cre- 
ation gap (i.e, the combined energy of a skyrmion- 
antiskyrmion pair) was found in the ideal 2D case 
to be exactly equal to half the gap in creating an 
electron-hole pair. (The latter is considered as an 
extreme case of a spin exciton with 5S Z = — 1.) 
Due to this result, at once the theory became by 
factor of 2 closer to the data found experimentally 
for the QHE activation gap. However, this and all 
later calculations, still remain in striking discrep- 
ancy with measurements. There is growing experi- 
mental evidenceE5T3 that at zero Lande factor, the 
energy required for activation of a dissipative current 
in a quantum Hall ferromagnet, is approximately a 
factor of 0.1 smaller than the calculated skyrmion- 
activation energy (the latter is one-half the creation 
gap). 

Justification for the application of the nonlinear 
0(3) model to the quantum Hall ferromagnet (QHF) 



is confirmed by a microscopical theory based on the 
Hartree-Fock (HF) Hamiltonian for the Coulomb in- 
teraction and on the approximation of-wave func- 
tions projected (WFP) onto a single LL.BB In these 
works, the energy of isolated skyrmionic excitations 
is recalculated and the minimum creation gap be- 
comes the same as in Ref. M Another approach, de- 
veloped by Iordanskii et alUdoes not use the approx- 
imation of WFP. The authors describe a skyrmion 
excitation as a smooth nonuniform rotation in real 
three-dimensional spin space. Due to the fact that 
the Coulomb Hamiltonian is invariant with respect 
to such a rotation, this approach has an evident 
advantage. The authors calculate the energy of 
skyrmionic excitations with the help of the pertur- 
bation theory technique. This theory uses for a bare 
Green function (GF) the appropriate mean-field one- 
electron GF. In doing so, the HF approximation js 
employed and the results (after a small correction! i 
turn out to be in agreement with earlier results. Qi3 
(See below in the Appendix II.) Naturally, any ef- 
fects of LL mixing are neglected there, and the re- 
sults reported in Refs. 0-p], as well as in Ref. El repre- 
sent the energy of skyrmions only in terms of a linear 
approximation of the Coulomb interaction; i.e. only 
in the framework of the first-order approximation in 
the parameter r c which is supposed to be small. 

In our work the energy of skyrmions is calcu- 
lated analytically with the help of the modified 
perturbation-theory technique. This technique is 
based on the excitonic representation (ER) which 
is suitable when a 2D EG is in a dielectric state, 
i.e. in the absence of free electrons and holes (see 
Refs. E9|-p6[). Neither the HF nor WFP approxima- 
tions are used. As in the case of Ref. M, a skyrmion 



excitation is considered as a rotation in 3D space: 



tp(r) = U(r)x(r), r = (x,y) . 



(1.1) 



Here -0 is a spinor given in the stationary coordinate 
system and x is a new spinor in the local coordi- 
nate system accompanying this rotation. The rota- 
tion matrix of U(t) (Ipty = 1) is parametrized 
by three Eulerian angles.EZI In the zero approxima- 
tion in terms of U(r) gradientsEa, we get x ex (J). 
Generally, in the limit of a vanishing Lande g factor, 
the results obtained may be presented as an exact 
expansion in terms of the parameter r c . Virtually 
two first terms of this expansion have been calculated 
exactly. 

As a first step, we ignore the Coulomb-interaction 
processes responsible for any decay of a one-exciton 
state due to the transformation into two-exciton 
states. Specifically, the part of the Coulomb Hamil- 
tonian kept involves: first, the terms responsible for 
the direct Coulomb interaction without any LL mix- 
ing; second, the terms providing the shift of the ex- 
change self-energy if an electron is transferred from 
one LL to another; and finally, the random-phase- 
approximation terms in which an exciton transforms 
into another one at a different point of conjugate 
space, but with the same spin state and the same 
cyclotron part of energy. In other words, such a re- 
duced Hamiltomian corresponds to a proper mean- 
field approach formulated in relation to excitons, but 
not to a quasiparticle excitation (as it would be in 
the case of the HF approximation). This Hamil- 
tonian is "excitonically diagonalized" (ED) and all 
one-exciton states (spin waves, magnetoplasmons 
and spin-flip magnetoplasmons V-aresent a full basic 
set for its exact diagonalizationBa The main idea of 
the present work is to employ a basic set of one- 
exciton states for the perturbation-theory formal- 
ism, instead of the basic set of one-particle states. 
If we restrict our study to the terms of the excitoni- 
cally diagonalized Hamiltonian (EDH), then the en- 
ergies of skyrmion excitations may be found in a very 
simple way for any magnitude of r c . In the QHF, at 
v = 1, and in the strict 2D limit we obtain the energy 
of the hole-like skyrmion E ED+ = jy / n/2r c huj c , 
and the energy of the electron-like anti-skyrmion 

is the skyrmion 
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(1.2) 



In the EDH model the gap thereby is determined by 
the smallest value among the cyclotron and Coulomb 
energies. At the same time, the EDH approach gives 
exact results of the first order in the expansion in 
terms of r c . In the limit of r c — > they are curi- 
ously in agreement with the resulte,obtained within 
the HF and WFP approximationsatl when the gap is 
i yjir /2r c fiuj c . This coincidence of the r c — ► results 



seems to be connected only with a special symme- 
try of the system studied (see the discussion in Ap- 
pendix II). 

The rest of the terms of the total Coulomb Hamil- 
tonian are responsible for other various LL mixing 
processes and provide additional corrections of the 
second and higher degrees of r c . The calculation 
carried out in Sec. IV (at v — 1) yields the exact 
second-order corrections E, ' = —O.OOSr^huic for the 
skyrmion and AE^ = ~0.382r 2 ,hui c for the creation 
gap. These corrections are independent of the mag- 
netic field. 

It is interesting that in the opposite case, when 
r c 3> 1, the Coulomb terms not involved in the 
EDH give again only small corrections of the order 
of hui c /r c to the fermion creation gap (see the end 
of Sec. IV). Therefore, if we formally consider the 
r c — > oo limit (by keeping the magnetic field con- 
stant we can study the e — ► limit), then for an 
ideally clean sample the EDH formula (1.2) repre- 
sents also the correct result A ED — ► Tiu> c . This has 
to take place for skyrmions which are characterized 
by a smooth spatial function on the length scale of 
Ib- The reasoning explaining this outcome may be 
as follows. Indeed, at a fixed v = 1 the parameter 
r c has one more meaning: it is the average inter- 
particle separation in units of the effective Bohr ra- 
dius as — h 2 £/m* e e 2 . When r c ^> 1, the average 
Coulomb interaction is weak as compared to the ef- 
fective Bohr atomic energy Eb = m*e 4 /H 2 e 2 , and 
it ceases being responsible for the gap determined 
by lowest-energy spatial excitations-tiJ However, in 
accordance with the Kohn theorem, Ell the cyclotron 
frequency remains always a relevant parameter of 
the clean system. If r c ^> 1, then Kuj c is the smallest 
quantity in the energy scale, hto c <C e 2 /eIb <C Eb, 
and the fermionic gap has to approach htv c or zero. 
Since the v = 1 condition seems always to provide 
an insulator phase of the system studied, the result 
A — > hto c is natural in the considered limit. 

Of course, these simple speculations ignore a dis- 
order which turns out to be the main factor, and 
it really governs the spectrum of the system just at 
r c 3> 1. This case is practically realized at com- 
paratively low magnetic field. Then the picture at 
large r c is determined by the competition between 
disorder and magnetic field or between disorder and 
interaction, and it turns out to be rather diverse (e.g. 
see Refs. p2| , p3| and the works cited therein). 

Meanwhile, in the clean limit just the situation 
when r c ~ 1 becomes experimentally relevant. Then 
our calculations within the EDH framework (Sec. 
Ill) as well as beyond of it (Sec. IV) demonstrate 
that the skyrmion-creation gap is substantially re- 
duced in comparison with the HF and WFP calcu- 
lations. In extrapolating to the region r c ~ 1, we 
can compare our results with the experimental data 
(Sec. V), if only for the highest magnetic fields at- 
tainable. 

We note that a mean field study of the r c ~ 1 fer- 
romasjiet was carried out numerically in the recent 
work.E-3 The statement of the problem used there 
seems to correspond to the EDH model, and the 
dependences found of the skyrmion energy on r c re- 



veal qualitatively the same trend as in the present 
paper. Unfortunately, a direct quantitative compar- 
ison with our results is impossible because in Ref. |34| 
the authors report only the results for finite 2D gas 
thickness or finite skyrmion spin. 

To conclude the Introduction, we comment shortly 
on the ER method used in the present paper. The 
ER technique implies a change-over from the Fermi 
creation operators, which generate one-electron 
eigen states of an ideal electron gas, to new exciton 
operators generating one-exciton states in the 2D 
electron system. When acting on the ground state, 
these exciton operators produce a basic set of exci- 
tonic eigen states. An essential part of the Coulomb 
interaction Hamiltonian (precisely the EDH) may be 
diagonalized in this basis. We extend in this work 
the range of the ER method and consider the basic 
set of two-exciton states. By using a special commu- 
tation algebra of the exciton operators EtEj (see also 
Sec. IV. A of the present paper), ER provides a sim- 
ple way to calculate matrix elements corresponding 
to various types of interactions. These may be not 
only the Coulomb terms (which are not involved in 
the EDH) but also e.g. electron-phonon or electron- 
impurity interactions. In terms of the ER, rthey 
are renormaliz cd . re&p ectively into inter-excitonEj'O, 
exciton-phononr^r 5 !! 26 ! or exciton-impurity interac- 
tions. In some particular cases the ER operators 
were used even in Ref. 09 when studying a two- 
component Fermi system with the symmetric model 
of interaction. Then it was found that this model 
corresponds to "inter- valley" waves in the 2D, semi- 
conductor at v = 1 in a high magnetic-field c3 or to 
spin waves under the same conditions .E3Ea 

II. VARIATIONAL PRINCIPLE 

We follow the general variational principle of E = 
minj /^\H\^\ //*|*\j . The averaging is car- 
ried out over the sample area. If we study an almost 
ferromagnetic state; i.e. the number of electrons A/" 
in the highest occupied LL differs from the num- 
ber Af<p of magnetic flux quanta by several units, 
| N — A/^l <1, then we can reformulate the varia- 
tional principle. The desired excitation presents a 
smooth non-uniform texture determined by the ro- 
tation matrix in Eq. (1.1). We divide the QHF area 
by the great number of Gi parts, which are much 
smaller than the total QHF area, but still remain 
much larger than the quantum of magnetic flux area 
2irlg . The energy of excitations of this type (includ- 
ing the ground state) may be found on the basis of 
the minimization procedure as follows: 



E = min 

u 



E 



mm 



(il>\H\tl>) 



Gi 



<v#>< 



(2.1) 



Here, the averaging is performed over a Gi area. All 
Gi areas add up to the total QHF area. The wave 
function tp should be substituted from Eq. (1.1). 

As to the outer minimization in Eq. (2.1), the 
only result required-ior its realization is the Belavin- 
Polyakov theorem.tZI Let us chose a unit vector n in 



the direction of the z' axis of the local system accom- 
panying the rotation. Evidently, n x = sin 9 cos <p, 
n y = sin 9 simp, and n z — cos 6, where tp and 9 are 
two first Eulerian angles. The minimum of th<f-$ra 
dient energy in the 0(3) non-linear er-model 



M 



min < — / d z r 



(d x n) 2 + (d y ri) 2 



where the topological "charge" is 



q T = — / d z vn ■ (d x n) x (d y n) 



47r|g T |, 

(2.2) 

(2.3) 



This corresponds to the degree of mapping of the 
2D plane onto a unit sphere of n directions and 
therefore, it is equal to integer number: q T = 
0, ±1, ±2,.... 

The procedure of the inner minimization in Eq. 
(2.1) is equivalent to the solution of the Schrodinger 
equation within the area Ax Ay = Gi, where the 
rotation is almost homogeneous, 

IdpUfAxAy « 1, \d 2 U\AxAy « 1, n = x,y. 

(2.4) 
At the same time, in zero approximation in terms 
of gradients, the "regional" state is a ferromagnetic 
with a great number of electrons corresponding to a 
local magnetic flux number 

A^(r) = AxAy/2nl 2 B > 1. 

For every region Gi the first and second derivatives 
dfJJ and d 2 u ll should be considered as external pa- 
rameters which depend only on the position of Gi 
(e.g. r is the position of the center of Gi). 

The substitution (1.1) into the Hamiltonian is re- 
duced to a trivial replacement of tp with \ in its 
Coulombjpart, but the one-electron part becomes of 
the form:LJ 



Hi 



24/ d2rf+(r) (" 



•iV 



f^) 2 x(r) 



(h = e/c=l, l = x,y,z), (2.5) 

where cr x ,y,z stands for Pauli matrices, and the 
parameters, il x ' y ' z (r) are proportional to small 
gradients CJ 

n* = i (l + cos*)^, 



£2„ = - (sin cos ifd^ip — sin ipd^O) , 



fi y = — (sin 9 sin <p d^ip + cos tpd^t 



(2.6) 



If we were to restrict our study to the one-particle 
approximation and neglect the Zeeman coupling, 
then the additional gauge field in Eq. (2.5) will not 
give any corrections to the one-electron energy. (In- 
deed, in this case one can turn every spin in any 
way without any change of energy.) However, this 



field changes effectively the "compactness" of the 
one-electron state at a certain LL, because an ad- 
ditional "magnetic field" V x (Q l aA appears. For 

electrons belonging to the upper occupied LL (at the 
filling v = 2m + 1 it has the index to) this additional 
field is 



SB = (x 



V x 



($01) 



x ) = v x n z 

Gi 



(2.7) 



and using Eqs. (2.6) we find also that it is equivalent 
to 

5B = -]-n- (d x n) x (d y n) . (2.8) 

The number of states within a LL is determined ex- 
actly in terms of one-electron wave functions. This 
value is changed by SNq — AxAySB /2n for the 
level to and the total number of states is changed by 
J d 2 r5B/2ir. Finally, due to the principle of maxi- 
mum filling, the topological injWiwt (2.3) takes on a 
new meaning microscopically rloBCl q T is the number 
of deficient (q T > 0) or excessive (q T < 0) electrons, 
i.e., 



TV = Aft - q T , where Af$ = ^~] N<f,(i 



(2.9) 



III. MODEL OF THE "EXCITONICALLY 

DIAGONALIZED" COULOMB 

HAMILTONIAN AND THE FIRST ORDER 

APPROXIMATION IN r c FOR 

THE FILLING FACTOR V = 2m + 1 

The additional field in the Hamiltonian (2.5) de- 
termines a certain perturbation operator Vh , and we 
can present the full Hamiltonian of the Gi region in 
the following form: 



H G = fl i + flint + Vh , 



(3.1) 



fixed g T .Thus even for results associated with the in- 
teraction part H mt of the unperturbed Hamiltonian, 
we have to take into account that the magnetic field 
is changed effectively by the value (2.7) and (2.8) for 
electrons within the LL to, and therefore the effec- 
tive magnetic length for this level is 



I, 



l B -l 3 B V x Q, z /2. 



(3.4) 



At zero approximation in terms of Vh, the ground 

state of this Gi region presents itself as the QHF 

with a total spin ^N<j, aligned along the z' axis of 

~ 2 
the local system, where N$ = AxAi//2ttIb ■ The 

Coulomb interaction does not change the spin of the 

ground state. If writing 



flint = fli 



ED 



ft 



int : 



(3.5) 



we extract from the Coulomb Hamiltonian the well- 
studied ED part (e.g., see Refs. |29] , |37| , and the next 
Section of the present paper) and will ignore so far 
the flint terms. The ground state of the Hamilto- 
nian floi + Aed is the same state |0) as it is for 
r c = 0. The ground state energy Eo m is determined 
exactly. In the case of v = 1 this energy is pro- 
portional to N,p T e 2 /sIb- Therefore the appropriate 
correction (within the ED approximation) is a 



SEoo = ^E ao l 2 B V x & 



for v = 1 . 



(3.6) 



(Here floo is considered after subtraction of the pos- 
itive background energy.) 

In the case of v — 2m + 1 > 3 the appropriate 
analysis reveals that we have to take into account the 
correction (3.4) only in Aed terms associated with 
the single level m, and also in the terms responsible 
for the exchange interaction between electrons of the 
m-th level and electrons of other filled levels having 
the same spin state. The result is 



where 



Vh 



2to* 



d 2 rj? 



over Gi 



-iV + A + Q l ai 



HV + A) 2 x, (3.2) 



and 



flint =|| rf 2 rd 2 r'x+(r)f(r') 



over Gi 

■U(r~r')x(r')x(v), (3.3) 

However, by using Vh for a perturbation technique, 
we should be accurate in avoiding a situation where 
we would be solving the Srodinger equation with dif- 
ferent numbers of electrons for perturbed and unper- 
turbed parts of the Hamiltonian. (Indeed, the num- 
ber of electrons depends on q T and therefore, on the 
perturbation term.) We will solve the problem at a 



6E 071 



-N^{e 2 /e)l B V x n z I m , 



(3.7) 



where 



J-m — 



rf 2 q 

(2tt)2 



V(q)e-i 2 / 2 P m (q 2 /2) . (3.8a) 



Here V(q) is the dimensionless 2D Fourier compo- 
nent of the Coulomb potential [in the ideal 2D case 
we have V(q) = 2ir/q, (here and everywhere below q 
is measured in units of l/Zg)], and P m (z) is a poly- 
nomial of the 2to power. If we set Po = 1/2, then 
the formulae (3.7) and (3.8a) determine the energy 
floo and the correction (3.6). For to > 1, we get 



Pm — r [im( z )J " 



A;' 

— -.z m ~ K \L 



-k 



-\z)\ 



k=o 



= - [L m {z)f 


+ i+X> 


t (z) [L k (z) - 


- Lk- 


-i(*)] , 






fc=l 






(3.8b) 



where L 3 k is a generalized Laguerre polynomial. [The 
simple derivation for Eqs. (3.7)-(3.8a,b) may be car- 
ried out in terms of ER by means of Eqs. (4.12)- 
(4.14) presented below the diagonal part of the Heu 
Hamiltonian.] 

To find the perturbation term Vq we substitute 
the expansion x = J2ap c ap<Pap into Eq. (2.5), where 
we chose the Landau-gauge functions as a basic set 
of functions cf> ap . The subscript p distinguishes the 
different members of the degenerate set of states, 
and the label a is a binary index 



a = (n a: a a ) , 



(3.9a) 



which represents both the LL index and spin index. 
Another designation will also be used when n or n 
is exploited as a sublevel index. In such a situation 
it means that 



n = (n, 1"), and n = (n, J.) . 



(3.9b) 



By integrating over the Gt area in Eq. (3.2) we 
should substitute x + £i and y + £2 for r-components 
and fl l + £1 d x il +£2^^ for Q (r), and then perform 
integration over £1 and £2- After routine treatment 
we obtain 



Vn 



ujc (u + u + + w + w+) , (3.10) 



where 



I 2 



E(« 



A7+|vx^Wn+i 



N n 



-JV*] + l B n z _ (A+ - K+) , (3.11) 
and 



£TQ+ 1 + nlQ + — . 

t n 7b+-l ^ n n _[_i 

(3.12) 
The following notation is used: 



Af = A/" T + A"x = J2 f^" + ^ 



where 



#n = E« 



^=E 



% p c„p ; 



^ = V Vn + lei— 65* 



'np ? 



n,p p 

(3.13) 
We employ here the designation a p (b p , c p ,...) for 
the electron-annihilation operator corresponding to 
sublevel a, (6, c,...). 



The sign of approximate equality in Eq. (3.10) 
means that we have omitted in the expression for 
U terms of the form F + a + + F-tT_ [where a± = 
(a x ± ia y )/2, the factors F± are of the order of 
(!bV) 2 ]. These terms are responsible for the devia- 
tion of all spins as a unit about the z'-direction, and 
they do not result in any contribution to the energy 
in absence of the Zeeman coupling.Ej 

The second sum in Eq. (3.11) corresponds to the 
formal change of the cyclotron energy due to the 
renormalization (2.7). Both operators N and K + = 
K^+K^ commute with the interaction Hamiltonian 
(3.3). (This feature of K + is a corollary of the Kohn 
theorenO; and in addition we have 

u> c K?\ •) In case |0) is the exact QHF ground state, 
then (Hoi + -Hint ) |0) = -EolO); and one finds also 

that (A^-AM |0> =JV^|0). 

First, we consider the correction determined by 
the U terms in Eq. (3.10). Suppose that we have the 
v = 1 filling. In this case A^O) = A+|0> = 0. The 

correction determined by U operators in Eq. (3.10) 
can be written in any order of r c in the general form 



Kl v H 0l 



SEi 



u c l B 



U 3 )' 



v x n z 



<0|AA|0) 



J 2 B V x?t z J2 n(0|iV n |0) 



n=l 



+ (coj B ) 2 N tP n z + n: 



-|A+|0) 



En-E+ 



(3.14) 



Here 



= Kf\0) is the eigen state of the total 
unperturbed Hamiltonian (the corresponding energy 
is E + = Eq + uj c ). Evidently, the first and third 
terms in Eq. (3.14) will always give a correction 
independent of r c . The second term also does not 
result in any correction of the first order in r c . [This 
sum gives only corrections of higher powers of r c 
which appear due to terms 7ii n t in the Hamiltonian 
(3.5).] 

The desired correction proportional to r c is 
thereby determined only by the operators W and 
W + in the perturbation (3.10). At the same time, 
by operating on the ground state |0), the terms H. m t, 
as well as the terms of the operator Hint X Vq , raise 
the cyclotron energy. Hence, the procedure of the 
perturbation theory in terms of 7ii n t would give only 
second- or higher-order contributions to the energy 
in terms of r c . 

Now, let us consider v = 2m + 1 > 1. By a sim- 
ilar analysis we can see that the operator (3.11) re- 
sults in a contribution independent of the Coulomb 
interaction or leads to other corrections which are 
of the order of r 2 and of higher orders in terms of 
r c . These corrections appear only on account of the 
Hint terms. If we restrict our study to the EDH 
model, then |0) = |0), and each of the operators 
Kf and Kt commutes by itself with Hed- These 

operators create the degenerate state Kf, |0) with 
energy Eq + lo c . Thus, if we want to solve the prob- 



lem to the first order in r c , and/or remain within 
the frameworks of the EDH model, then we may use 
the U terms only to obtain the zeroth order contri- 
bution to the final result. (Such contributions from 
all terms of Vq cancel each other in the zeroth order 
of r c .) 

In this section, we consider the EDH as a Hamil- 
tonian responsible for the Coulomb interaction. As 
we have seen the operator (3.12) should really be 
taken into account as a perturbation. Moreover, 
only operators Q + — with n = m— 1 and n = m 
have non- vanishing results of operation on |0) after 
their commutation with Hed- Therefore, only the 
last term of the operator W contributes to the en- 
ergy of skyrmions. 

A. Filling factor v = 1 (m = 0) 

In this case, the state 



Therefore at |g T | = 1 the creation gap is equal to 



(1 + Wc) 



(3.20) 



(this is the factor before |q T | /2). Within the strict 
2D limit, when V(q) — ► 27r/g, we get 



CSF 



el B ~" 2el B V 2 ' 
and we arrive then at the result in Eq. (1.2). 

B. Filling factor v = Ira + 1 (m > 1) 

In this case, there are two basis states Qi_ 1W r |0) 

and Q + — -JO) forming spin-flip magnetoplasma 

modes (with zero wave vector) of the unperturbed 
Hamiltonian. In the 2x2 matrix equation of 



|sf) = q+ t 



(3.15) 



is an eigen state of the unperturbed Hamiltonian 
which corresponds-itei the so-called "spin- flip magne- 
toplasma" mode E3e3 with a zero wave vector: 



Hoi + H ED: Q+ T |0) = (oj c + 8 s l) |SF) . 



Here, £q± is the Coulomb part of the energy: 



£ s l = (e 2 /d B ) 



d 2 q 



q 2 V(q)e- q2/2 (3.16) 



If one sets formally £qi = 0j then the first- and 
second-order corrections of the one-electron energy, 
in terms of the perturbation Vh, are exactly can- 
celed in the result. [One can check this fact with 
the helpof Eq. (3.14) and by employing the useful 
identity!] of V x & = 2eyfeO : ' x fi fe .] Thus, we obtain 
the second-order correction determined by the w c W 
perturbation: 



8E- 



2 Ax Ay + 

W 2 = ^c—^ ^I^ + 

Z7T 



(sf| q+ t 



i 



U! c + £q1 



(3.17) 



The factor 17_r2 + can be expressed in the terms of 
the unit vector n, since 



Hed,Q 



jjfi 



k=m — 1 



(j = m-l,m), 

the diagonal elements are 



(3.21) 



C 33 ~ £ J 



rf 2 q g 2 ^(g) c - g V2 r £ i ( 2 /2)1 - 
(27r) 2 2(j + l) L^W/^J 



(j = m-l,m). 



(3.22) 



The off-diagonal matrix elements are equal to each 
other: 



c(0) _ _ c(0) 

m-1 m m m+ i 



d 2 q g 2 V(«) 



(2tt) 2 2 1 /m(m+l) 



xe^ 2 / 2 i r 1 „_ 1 (< Z 2 /2)i,i l (g 2 /2). 

Two spin-flip modes thereby have states with ener- 
gies co c + £ ± F , where 



£T = 



eh 



(Cm + £m-l) /2 



±Y (e m - e m -i) 2 /4 + w 2 



(3.23) 



0+0+ = i 

T ± a 



(^n) + (<9 y rf) =F2n • (c^n x <9 y n) 



(3.18) 
After integration over the 2D space (<i 2 r = Air Ay) 
we obtain with help of Eqs. (2.1)-(2.3) and (3.6)- 
(3.8a,b) the skyrmion energy corresponding to the 
charge q T : 



3f 

2el, 



-loq-i 



(3.19) 



If we were to neglect the values of the commutators 
(3.21), then again all corrections determined by the 
perturbation Vh will add up to zero. The non-zero 
result for the Gi region is determined by the second- 
order correction of the perturbation theory and is 
caused by the W operators 



SE W 2 — 



AxAy 0+ 

II 11 , uj r r c 

2tt " + c 



(m + l)e m + memr-i + r c (2m + 1) (e m e m _i - w^) 



l + r c (e m + e m -i) + r% (e m e m _i - w^) 



(3.24) 



After summation of the combination SEyy2 r +SEo m 
over all such Gi regions, we find the energy of the 
skyrmion. Using Eqs. (2.1)-(2.3) and (3.18) we see 
that in the first order in r c it takes the form 



-Eed(<7t) — —r- 
elB 



e 2 f (I«t| — g T ) 



[(m + l)e m + me m _i] 



3 r 



(3.25) 



Therefore, in this r c — > limit, even for m = 1, 
the skyrmion-antiskyrmion creation gap is essen- 
tially larger than the gap for the electron-hole pair. 
Indeed, with help of Eqs. (3.8a,b) and (3.22), in the 
ideal 2D case we can obtain the creation gap which 
turns out to be equal to 

_e^_ll R 

d^"8~V 2 ' 

whereas the appropriate value for the electron-hole 
pair is 

d B 4 

Just the latter determines thereby the activation 
charge gap in QHF at v = 1. Analogously, one can 
prove that for any v > 3 the skyrmion gap found 
from Eq. (3.25) is larger than the quasiparticle gap. 
Thus, skyrmions are lowest-energy fermionic exci- 
tations only in the case of the filling factor v = 1 . 

IV. CORRECTIONS AT V = 1 TO SECOND 
ORDER IN r c 

Generally, the second-order correction to the EDH 
skyrmion energy E ED is a combination of several 
parts which have different origins. 

First of all, when calculating the ground-state en- 
ergy to zero order in Vh (but to second order in 
r c ), we must again take into account the renormal- 
ization (2.7) and (3.4). The corresponding value, 
being of the order of A^r 2 w c , is negative (as it has 
to be for any second-order correction to the ground- 
state energy). Therefore, the renormalization cor- 
rection turns out negative at positive SB, namely: 
SE 0r 2 ~ — A^V x Ct z r^.u! c . Precisely the same form 
of correction we obtain for SEy^ which is caused 
by the second term in Eq. (3.14). However, this cor- 
rection is surely positive in the case of V x il z > 0. 
In the following, we will see that both corrections 
cancel each other: 



SE, 



Or; 



SE. 



UH 2 



0. 



(4.1) 



Another correction of the required order is deter- 
mined by the fourth order of the perturbation theory 

in terms of the sum ii- m t + <*> c ( VV + VV ) . This cor- 
rection is quadratic in 7ii nt and in VV and should 
take the form 



SE 



W 2 H 



2 = N<j,r c oj c 



( m n1:n+ + m nzn+) (4.2) 



(l?7i| ~ I V2 1 ~ 1)- For the calculation of the cor- 
rections studied the perturbation technique can be 
formulated in terms of the ER. 

A. Excitonic representation 

We proceed from the following form of the inter- 
action Hamiltonian (cf. Ref. 23): 



Hint = -7T- 2_^ V M ca(q) 



a,b,c,d 



x exp [iq x (p' - p)}b+ +qy d+c p , +qy a p , (4.3) 



where 



Hdca(q) 



e 2 V(q) 



h n b n a ( C i)K c n d ( f l)^ a ,a b Sa a ^ d ■ 

(4.4) 



2tteIb 
The function /i„fc(q) is 

where 



q± =T-^(q x ±iq y ) 



(4.5) 



In the ER we change from electron cr ea ^ iq n (anni- 
hilation) operators to the exciton ones:li^ 



Qabq = Qt a _q W 



(4.6) 



This is a generalization of operators (3.13) in the 
case of non-zero wave vector q. A one-exciton state 
is defined as 



Kq) = Q: bq |o). 



(4.7) 



We will also use the intra-LL "displacement" oper- 
ators 



A = N^/ 2 Q+ , 



(4.8) 



for which, evidently the following identity takes 
place: 

A+|0) 

<5o,q|0) if n a < to, or n a — m and cr a = +1/2, 
if n a > to, or n a — to and <J a — —1/2. 

(4.9) 
The commutation rules of the operators (4.6) and 
(4.8) present a special Lie algebra (cf. Refs. |2l] |23| ): 



-^abqi i -^abc{2 



W-a6qi j W,a6q2 



2fccqi ' 2a 



b(l2 






= o 



(a^b^c^d), (4.10a) 



^6cqi ' ~^abci2 



= V /V * (qiXq2W2 2ac qi+q2 (c ^ a, 6), (4-lOb) 



Wlabqi ; ^ahq 2 



_ e *(qixq2)z/2^ 



qi-q 2 



-i(q 1 Xq 2 ) i /2 



5, 



qi-q 2 i 



(4.10c) 



i(qiXq 2 )z/2 4 Q _ -i(qiXq 2 )z/2 n £) 



ihq 2 



= AT- 1 Q a6qi+q2 . (4.10d) 

The interaction Hamiltonian (4.3) may be rewrit- 
ten in the form 



q,a,b,c,<i 



^ ^aba(q)S+ 



(4.11) 



[In particular, if a = (j, j) and b = (j + 1,4), then 
£ab(0) = £J°\ see above Eq. (3.22).] 

However, if to > and Sn > 0, then the EDH also 
involves an off-diagonal part. When operating on 
the state (4.7), the off-diagonal terms give a finite 
combination of other one-exciton states. Thus, 



H 



ED 



H ED 



ErVoff-di 
H ab 

ill.) 



(4.16) 



Contrary to the definition (4.12), the summation in 
Eq. (4.16) is carried out only over the ab pairs in 
which the sublevel a is occupied and the sublcvel 
b is empty in the state |0). The members of this 
summation are 



rroff-di 
H ab 



E 

l+n d =n b +n 



2 , \Vadcb{q)Q cdq QabcL 



VadbMQtavQdbd ■ (4.17) 



One can check that 

6q 



rroff-di Q+ 
n ab > ^ab 



|0>= E ^RVqlO), (4-18) 

a' b'^ab 



Now we can extract from this expression the ED 
part. At least these terms do not change the 
cyclotron energy. (In other words, they have to 
commute with the one-electron Hamiltonian Hoi.) 
Therefore, to find the EDH we should consider in 
Eq. (4.9) only terms with n a + n c = rib + n d- A part 
of these constitute an operator in which the states 
of the type of Eq. (4.7) are the eigen states. Such a 
diagonal part of the EDH can be written as 

#ed = E^ + E #°*> ( 4 - 12 ) 



M^»«<»l,) 



where 



H* = \ E U M (W*i+iq - ij) , (4.13) 



and 



H ab = ]T [Uab(q)N^A+B^ + U ab (q) (Q+ fcq Q abq 



and the pairs of the states a'b' ... in Eq. (4.18) pro- 
vide the same Sn and SS Z as those in the case of the 
pair ab: 



Sn = n b - n a = n v - n a < 



SS Z = a h — a a = g v 



(4.19) 



The finite set of equations (4.15) and (4.18) 
thereby determines the eigen energies and the eigen 
states of the EDH which correspond to given Sn, SS Z 
and q. [Specifically, in this way the spin-flip modes 
at q = (3.23) have been found.] 

All other terms of Hi nt = H mt — Hed, with which 
an operator Q^ b does not commute, have the fol- 
lowing form: 



H 'ab 



Z^ Z^ v s 

cd^af} q 
9^a 



Vgacd{q) ^dgq Qacq 



Y^ ^2Vcdbg{q)Qg C qQdbci- 



(4.20) 



,d^a,b q 



-B- 



(4.14) 



We have used the notations U a b = V a bba an d U a b = 
Vbaba- One can check that for every operator (4.6) 
we get 



H a + Hb + H a b + 2_^ \Hac+ Hbc) , Q a fc q 
c^a,b 



|0) 



If operating on the state (4.7), these terms lead to 
"superfluous" two-exciton states. Meanwhile, some 
operators (4.20) do not change the cyclotron energy 
and even within the approximation of the first order 
in r c , they must be considered for the correct calcu- 
lation of exciton energy. Specifically, for the spin-flip 
mode (a = 0, b = 1, if v = 1) the terms 



^T = E F i°io(g)Qr!i q Qo 



Iq 



H. c. (4.21) 



CMQtJO) (4.15) 



also have to be taken into account as well as those 
of .Hed- 



We will calculate the second-order corrections to 
the energy in the case of the filling factor v = 1. 
Therefore, within the framework of our problem, we 
are interested in results of the operation of TL ln t on 
the state |0) or <2~J- |0), and it should be chosen in 
the form 



Hi, 



E Vbd 00 (q) Qobq Q-Od -q 

q,nb>0,n d >0 



+ E Wq)2o + hq 2^ + H. c. (4.22) 

q.n b ,n d 

The terms (4.21) enter into the second sum of this 
expression. 

B. Two-exciton states at v ■=■ 1 

The operation of the Hamiltonian (4.22) on the 
EDH ground state |0), at v = 1, leads to two-exciton 
states of the type of 



\a)=P+\0), 



(4.23) 



where we will denote as P+ the two-exciton creation 
operator 



p+ = -0+ o + 

a c\ ~^Oaiq a -^Qa-z— q Q 



(4.24) 



and designate as a the composite index a = 
(ai,a 2 ,q Q ) [correspondingly j3 = (h, 6 2 ,q/3),-], 
which obeys the evident condition 



\a) = |ai,a 2 ,q Q ) = \a 2 ,ai,-q a ) 



(4.25) 



Any state (4.23) is a "quasi" eigenstate of the un- 
perturbed Hamiltonian H i + i?ED: since 

Hqi + H EU , Pa |0) = K (n ai + n a2 ) + £ 0ai (q a ) 



+S 0a2 (q a )}P+\0)+ \i) (4.26) 

els 

[£o ai and £oa 2 are Coulomb energies determined by 
Eqs. (4.15)], where the state \e) has a norm of the 
order of l/N^. However, in comparison with the 
set of orthogonal one-exciton states (4.7), the states 
(4.23) are "slightly" nonorthogonal to each other. 
We can find that 



^ai,6i ^a2 ,62 



J qa,q,3 



J(q a >np) z 



0ci2,bi 0a\ ,62 



J q Q ,-q,3 



i(q^Xq Q ) z 



N A 



(4.27) 



Meanwhile, this nonorthogonality has to be taken 
into account if we are to consider a combination 



Ewi/ 3 ) 




(4.28) 



(this is a summation over all of the members of the 
composite index). In this case, the function ipp — 



<p(bi,b2,q_p) turns out to be non-single-valued one. 
Indeed, let us project this state onto a certain state 
(4.23). We obtain 



where 



& 



Vn 



2 ( * Q 



¥>J 



by- 



PVP 



(4.29) 



(4.30) 



is a Fourier transform determined by the kernel 

T a& = <5 0lj6l <5 02 ,6 2 e^« xt ^* . (4.31) 

If <p a — Tp a , then any projection (4.29) is equal to 
zero. Only the "antisymmetrized" part 



(a) 



1 



[<Pp - Vp) 



(4.32) 



contributes thereby to the combination (4.28). The 
origin of this feature of states (4.28) is related 
to the permutation antisymmetry of the electron 
wave function in the system studied (cf. for ex- 
ample Ref. p9). Note also that Tp a = ip a and 
J2 a w(a 1 ,a 2 )^f a = Y, a w ( a i> a 2)^aWa~, where w 
is such that u>(ai,a 2 ) = tu(a 2 ,ai). As a result, we 
get the useful equivalence 



^w{a 1 ,a 2 W a 



*0) (a) 



^ W (a l7 a 2 )C^ a) 



^ W ( ai ,a 2 )C (a Va. (4.33) 



In terms of these definitions the expectation (4.27) 
may be rewritten as 



- ^ ( <W - Jf E Fa-yS-yrf J 



— °a,/3 ' 



where 



(4.34) 



^a,/3 — q \*ai,bida2,b2 < 3ci a ,ci/3 ~r Va 2 ,bidai,b 2 ^q ai — Q/3 J ■ 

(4.35) 
C. Perturbation-theory results 



When Hint +uj c VV + W + J is a perturbation and 

Hqi + Hed is an unperturbed Hamiltonian, then it 
is sufficient to employ as a basic set the two-exciton 
states (4.24) and the spin-flip state (3.15). Thus 
the correction to the EDH ground state |0) may be 
presented in the form 



S\O)=C o \0) + Y,Ca\a)+D\SF) 



(4.36) 



a^O 



where the factors C a and D, should be found in 
a specified order in terms of Hint and VV (actually 
in terms of r c and 0-). In our case, where we are 
concerned only with the antisymmetrized functions 
C a = Cq , the above equations (4.29), (4.33) and 



(4.34) show that two-exciton states may be considered as an orthogonal and normalized basis, for which the 
perturbation-theory technique may be used in its traditional form. 
First, we find C a to the first order in Hint, 

C ow = 0, C aH = -{a\ii- lnt \o)/^ c a (a^O). (4.37) 

where A^ = uj c (n ai + n a2 ) stands for the cyclotron part of energy of the two-exciton state |a) [c.f. Eq. 
(4.27)]. When substituting |0) = |0) + J2 a Ca«|a) into the second term of Eq. (3.14), we obtain 5E UH 2 = 
V x ffwjfl Sa ('Vl + n a 2 ) |Cqw| • Then by calculating also the second order correction SE 0r 2 = V x 
& z l% Xla ^aH (O Hint ct), we come indeed to the result of zero in the combination (4.1). 

The desired correction (4.2) is determined by means of a conventional procedure. In which, factors C a 
have to be found sequentially up to the second order in VV and to the first order in Hint- Whereas D, is 
determined to the second order in Hint and the first order in W. The result is written in the form 



SE W 2 H 2 =Ai+ 2A 2 + A3 



where 



(4.38) 



(4.39a) 



(4.39b) 



(4.39c) 



(we set here A| F = uj c which is the cyclotron part of energy in the state |SF)). The matrix elements entering 
into these expressions are calculated in Appendix I. 

Consider for example, the term A\. With help of Eqs. (Al.l), (Al.2) (A1.5), and (A1.7) in the Appendix 
I and using Eq. (4.33), we find that Eq. (4.39a) is changed into 
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A\ = -r 2 c u c ~N <i , 



n+n- 



9 a 9a {n a 



Z^ 



(n ai + n a2 ) (n ai + n a2 +l) 



Q+n- 






* ( a ) 1 1 -v 

ffaffa ("«i +n a2 ) 

{n ai + n a2 ) (n ai + n a2 - 1) 



(4.40) 



After substituting Eq. (Al.2) for g a , the suitable sequence of mathematical treatments is as follows: we 



perform the summation over all of n ai > 1 and n a2 > 1 keeping the sum n a = n ai 



> 2 fixed; then we 



make the integration over q a [the antisymmctrized function <?q already contains an integration according 

to Eqs. (4.30) and (4.31); therefore the second term in the expression g^Qa — (|<?a| 2 — 9a9 a )/^ leads to 
twofold integration over q,3 and q Q which, however, can be reduced analytically to a simple onefold integral]; 
and finally the numerical summation over n a is performed. In the ideal 2D case [i.e. if V(q) = 2ir/q] the 
result is 



Ai = -rluJcN^ (0.0056fij:^ 



0.0077ftjftl 



(4.41) 



In a like manner, the calculation of A 2 and A can be carried out. In so doing, for the ideal 2D system 
limit, we obtain A2 = — O.O248r 2 ,w c A0$7+Sl+ and A3 = — 0.318r^.u} c N^iT^i}^.. [Notice that the operators 
(4.21), which do not change the cyclotron energy, contribute only to the term As-} 

With help of Eq. (3.18) and Eqs. (2.1)-(2.3) the summation over all of the Gi regions yields the second- 
order correction to the EDH result (3.19): 



E 



(2) 



= E 5E w 2 n* = r\ 



y {m-m) + — (vi 



m) 



(4.42) 



where 771 = —0.374 and 772 = —0.0077. The correction to the EDH skyrmion-antiskyrmion creation gap is 

AEW =r%w c {r h + Vi) • (4-43) 

With Eqs. (3.19) and (4.42) we find the total r 2 , correction to the HFoEl result: 



E, 



(2) = 

tot q T 



(Ml 

2lOr 



(2) 



(\q T \-q T )+E%. 



The corrections for quasiparticles are correspond- lating the correction to second order. Needless to 



ingly: E i 



(2) 



-0.767r 2 u> c for the electron-like anti- sa Y ii; would be incorrect to apply this r 2 c correction 



AE t 



(2) _ 



tot + 



-0.0077r2 Wc 



for the holclikc 



skyrmion 
skyrmion 

To conclude this section we prove that all per- 
turbative (in terms of Hint) corrections to the EDH 
gap (1.2) vanish in the r c — v oo limit. Note that 
the above calculations of SE 0r 2 , SEyj-p and EVv 2 w 2 
present the corrections of the second order in Hint 
where the energies of one- and two-exciton basis 
states are considered within the zero approximation 
in r c . At the same time, the technique used pro- 
vides a formal possibility to develop perturbatively 
an expansion in Hint for arbitrary, r c . To do this we 
should replace Asf and A^ - in Eqs. (4.37) and 
Eqs. (4.39a-c) with their values exactly calculated 
within the EDH model. That is we should add the 
corresponding Coulomb shifts: Asf — * ^c + ^oi an< ^ 
A c a -> uj c {n ai +n a2 ) + £Q ai (q a ) + £oa 2 (q a ) [A^ 7 = 
A c a (a^b,c), see Eqs. (3.16) and (4.15)]. Naturally, 
at r c < 1 this procedure becomes senseless if we were 
to restrict our consideration to the second order in 
Hint only. Indeed, both of the operators Hed and 
Hint are proportional to e 2 /eIb and the accounted 
Coulomb shifts in the denominators of Eqs. (4.37) 
and (4.39a-c) would only yield the r% and higher 
order corrections which are beyond this approxima- 
tion. 

However, let r c > 1, and estimate at once all terms 
of the perturbative expansion in H n t ■ In this case 
\C a \ ~ 1 to any order in H n t, and now there is 
no cancellation (4.1), because SEq t ^ ^ —SEuh* 
for any k > 1. (The more specific estimations are 



5E 0r K ~ VxCl z l%u! c r c and 5Euw 



Vxtt z l 2 B uj c 



) 



Nevertheless after integration over 2D space both of 
these corrections contribute only to the term pro- 
portional to the charge q T (2.3), and therefore they 
do not contribute to the creation gap of skyrmion- 
antyskyrmion pares. As to corrections SEy^2- H K , we 
find easily that they are all of the order of l B ui c V 2 /r c 
and give a correction to the gap of the order of 



K>2 



AE™ ~ uj c /r c 



When r c ^> 1, the EDH value Aed ~ oj c presents 
thereby the main part of the creation gap. 



to the case when rj, > 1. Nevertheless, the perturba- 
tion theory result indicates at least the tendency of 
the creation-gap variation with r c . Fig. 1 displays 
also the corrected value 



A- r = A E 



A£< 2 



where 



AE^ = ( Vl + m )r 2 c tu c 



(5.1) 



(5.2) 



We choose conventionally the region r c < 0.2 as a 
region of the small r c values where the result (5.1) 
is correct. The curve for A^™ r trends thereby to a 
more severe decrease in the gap. In addition, we ex- 
tend this curve further (by the dashed line) approxi- 
mately to x« which corresponds to the experimental 
conditionaij ( the creation gap measurements made 
in Ref. [13] seem to be the most high magnetic field 
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FIG. 1. The normalized creation gap of the ED Hamil- 
tonian model and the corrected gap are shown above. 
The dotted line represents the HF approximation result. 
The experimental data (closed symbols) are from Refs. 
3 (the square), 4 (the circle), and 5 (the triangle). 



V. DISCUSSION 

Thus, our calculations consist of two main stages. 
In the first stage we have considered only the ED 
part of the Hamiltonian, where the LL mixing is 
partly taken into account. The corresponding cre- 
ation gap for charged quasiparticles A BD (i.e. for 
|<7 T | = 1 skyrmion- antiskyrmion pairs) is shown in 
Fig. 1. Here, r c is an arbitrary parameter (formally 
it does not need to be small in this approach). We 
see that even the EDH model reflects a significant 
reduction of the gap with a growing r c . Besides the 
obtained A ED yields the correct limiting values for 
r c — > and for r c — > oo. 

In the second stage we have treated the remaining 
part of the Hamiltonian perturbatively in r G , calcu- 



Such results for A ED and for A^ r demonstrate a 
significant role of the LL mixing. The curves in Fig. 
1 become closer to experimental data (c.f. the HF 
value shown by the dotted line). We demonstrate 
here only the data found at fairly strong magnetic 
fields, namely: B = 11.6, 5.5 and 3.5 T corre- 
spondingly in Refs. |l^Jl^,|l^. Meanwhile, it should 
be noted, that the creation gap experimentally ob- 
served depends sharply on the effective g factor in 
the vicinity of g = [in accordance with the low 
(glng) 1 / 3 , see Ref. [jj. The result presented in the 
work of Ref. |l3J and marked here by the square is 
really obtained at zero Zeeman energy (tp_.aii ex- 
perimental accuracy). The measurementsE3il3 are 
carried out at small, but still not vanishing g fac- 
tors. Therefore, we have extrapolated these data to 



the point g = 0. The results of this conventional 
extrapolation are presented in Fig. 1 by the closed 
triangle and the circle . 

At the same time, the dashed line for value A^ r 
evidently presents an underestimated result. In the 
higher order corrections in terms of r c the curve cor- 
responding to the true energy gap should pass some- 
where between the curves of Fig. 1 (the smaller the 
parameter r c is, the closer this true gap should be 
to the calculated value A^™ 1 ). In the ideal 2D case 
the results A BD and A"£™ calculated for r c ~ 1 may 
be considered as upper and lower limits for a real 
creation gap. 

Finally we should emphasize that at least two im- 
portant effects have been ignored with the calcu- 
lation of the Fig. 1 curves: finite thickness (FT) 
correction, and disorder broadening of Landau lev- 
els (LL's). Both of these reduce the energy gap. The 
usual way to take into account the FT is to modify 
the Coulomb interaction: V(q) = 2irF(q)/q, where 
the formfjactor F is parameterized by an effective 
thickness. c3 Any second-order correction in terms of 
r c would involve this formfactor doubled and there- 
fore, it would be more sensitive to the FT effect. 
Roughly speaking: a 30% reduction due to the FT 
correction to the spin- flip- mode energy (3.16) causes 
the corresponding correction to A ED [see Eq. (3.20)] 
and determines a reduction by w 50% in AE^ [also 
by a factor of ~ (0.7) K in the following AE^ cor- 
rections] . The two curves A ED and A E " r would start 
thereby at r c = from m 0.7a/7t/8 and become more 
sloping. At any r c the gap turns out to be smaller 
because of the FT effect. 

The disorder may govern QHF features critically 
and even in the ground state this can lead to a 
realignment of spins with respect to one another. 
Specifically, the calculation of disorder effects de- 
pends on the model for the random potential. The 
white noise potential (arising e.g. due to chargeless 
point defects available in the 2D channel) is consid- 
ered in Ref. E2l A perceptible change of the charge 
gap seems to be related to appearance of a skyrmiopj 
like structure in the ground state. The authors^ 
found that the latter occurs at v = 1 starting from 
some appreciable threshold for the amplitude of the 
disorder potential correlator. In the opposite case of 
long range potential fluctuations (mostly determined 
by charged impurities situated out of the spacer), 
the gap should change smoothly with correlator am- 
plitude and could be estimated as follows. When 
chargeless exciton exists, the disorder broadening 
determines a finite cut-off value q m for 2D momenta: 
q < q-m- This momentum q m is related to a certain 
distance l^q m (in usual units). At this distance a 
force of the Coulomb interaction between quasipar- 
ticles that form the exciton, becomes equal to an 
external random force appearing due to the disor- 
der potential. Hence, the real creation gap for free 
quasiparticles decreases by a value E x (oo) — E x (q m ), 
where E x (q) is an appropriate exciton energy calcu- 
lated within the clean limit (c.f. the analysis in Refs. 
|23|Ji6|) . It is rather difficult to estimate the gap re- 
duction corresponding to our specific case, because 
the energy with q of a skyrmion-antiskyrmion exci- 



ton is unknown. However_£or the spin excitonHI 
the analogous estimations results in a reduction of 
« 20%, if the random force is caused by distant im- 
purities. (In real 2D structures this force could be 
estimated as ~ 0.1 K/nm.) 

Thus, the disorder and FT effects also play a role 
in the gap reduction at r c ~ 1. However, for the up- 
to-date 2D structures they seem to be less important 
compared to the basic effect of the LL mixing. 
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APPENDIX I: MATRIX ELEMENTS 

The commutation algebra (4.10a-d) for exciton 
operators (4.6) and (4.8) allows us with the help of 
the rule (4.9) to calculate the relevant matrix ele- 
ments with relative ease. Using Eqs. (4.22), (4.4), 
(4.5), (4.23)- (4.25), (4.29), (4.30), and in view of the 
fact that Q+ fi = Q^q we find 
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-LT 2 1 (^/2), if er 0l = L er a2 = -i. 



In the case of the state |/3) corresponding to 6S Z = 
(i.e. at,! = &b 2 = 1/2) and in view of Eqs. (3.12), 
(4.24), (4.34) and (4.35) we obtain 
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where 

d a/3 = fit (\/n bl + 15 a .j3 1+ + y/n b2 + l8 a ^ 2+ \ 



+0| (^ri^Sa.js^ + y/nb^8 at/ 3 2 _) , (A1.6) 

and indexes f3\± and /?2± designate the states 

\Pi±) = \nbt ±l,nb 3) qfi) , |/3a±) = |n 6l ,n& 2 ±1,^) 

[See definitions (4.35), (4.25), and (3.9).] In case 
n ai +n a2 = n bl +n b2 , we find also that for any function 

V a , 

^w 7 d* Q d 7 ^ = 5 a ,f3 [£L±Cl+vp 1+ (n bl + n b2 + 2) 
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APPENDIX II: 

OF THE TV — > 



THE COINCIDENCE 

RESULTS at 1/ = 1 



We investigate the origin of the coincidence 
of our results, which are exact at r c — + 0, 
with the results obtained:,, (i) within the HF 
and WFP armcoximations;Q (ii) within the HF 
approximatiorOQ. In any case the gap has to be pro- 
portional to e 2 /els- However, generally, the specific 
factor should be different in these three approaches. 

In the worku, where only a single LL is consid- 
ered, the corresponding factor is determined by the 
q 2 term of the expansion of £ sw (q) « q 2 £ s '( v (0)/2 at 
a small wave_v.ector q. Here S sw (q) is the energy of 
spin excitorSa. [This value is equal to £ o(g) in 
our notations, see (3.9b) and (4.15)]. The creation 
gap in the work of Ref. |5j turns out to be equal to 
the inverse spin-exciton mass M,," 1 = £^.(0). 

In the work of Ref. both n = and n = 1 LL's 
are used for the presentation of the bare one-electron 
GF. In so doing, the external HF field for this GF has 
been taken into account. Therefore, the denomina- 
tors of the GF contain the energies of one electron 
placed at the n = 1 LL in its spin "up" and spin 
"down" states. The result for the skyrmion-creation 
gap turns out to be proportional to the difference of 
these n = 1 energies. At the same time, the electron 
energies can be measured from the energy of a dis- 
tant hole at the n = level. Then, the required dif- 
ference corresponds to the extreme case (i.e. to the 
q — ► oo limit) of thc-difference of the corresponding 
cxciton energies OEJJ E x = [£ j(q) - £oi(q)] \ q _ loo - 

[The energies £oi(?) and ^ot(^) ma y be found di- 
rectly from Eq. (4.15).] 

We remind that our r c — v result for the gap is 
£pf = £ o y(0). Thus, generally, these three different 
approaches should lead to the different results. How- 
ever, due to specific features of the QHF studied, all 
these three values actually are equal to each other: 



M sw = E l - £ ll = -f- 
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oo exchange energies for the 01-magnetoplasmon. 
Therefore, the Eqs. (A2.1) appear to be nothing 
more than an coincidence peculiar to the system 
studied. 

If we study a single skyrmion or antiskyrmion, 
then we see that their energies (3.19) are determined 
also by an additional correction proportional to q T 
(where q T = ±1). In fact, this correction in the 
present work as well as in Refs. [?],|| is determined by 
the rotation-matrix feature (3.18) and by the renor- 
malization rule (3.4). Therefore, under the coinci- 
dence condition (A2.1), we arrive again at the iden- 
tical results. The approach of Ref. seems indi- 
rectly to contain also certain features analogous to 
(3.4) and (3.18). Also it results thereby in the same 
energies of an isolated skyrmion or antiskyrmion. 

Finally, it should be noted that for the filling v > 3 
our result (3.25) differs from the result of Ref. a 
This fact reflects the role of low lying LL's which 
participate in the skyrmion formation. Neverthe- 
less, the v > 3 skyrmion creation gap, just as in the 
approach adopted in Ref. pL turns out to be lower 
than the corresponding quasiparticle gap. (In the 
work of Ref. M only the v — 1 case was studied.) 



In particular, the coincidence of fgf and E\ is the re- 
sult of the "accidental" equality of the q = and q = 
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